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Abstract Helioseismology has been widely acclaimed as having been a great 
success: it appears to have answered nearly all the questions that we originally 
asked, some with unexpectedly high precision. We have learned how the sound 
speed and matter density vary throughout almost all of the solar interior - 
which not so very long ago was generally considered to be impossible - we 
have learned how the Sun rotates, and we have a beautiful picture, on a coffee 
cup, of the thermal stratification of a sunspot, and also an indication of the 
material flow around it. We have tried, with some success at times, to apply 
our findings to issues of broader relevance: the test of the General Theory of 
Relativity via planetary orbit precession (now almost forgotten because the issue 
has convincingly been closed, albeit no doubt temporarily), the solar neutrino 
problem, the manner of the transport of energy from the centre to the surface of 
the Sun, the mechanisms of angular-momentum redistribution, and the workings 
of the solar dynamo. The first two were of general interest to the broad scientific 
community beyond astronomy, and were, quite rightly, principally responsible 
for our acclaimed success; the others are still in a state of flux. 
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1. Prelude 



In the early heady days of helioseismology, the new techniques of inference, 
coupled with pertinent observations that had been stimulated by them, went 
hand-in-hand with consequent scientific discovery, and helioseismology was not 
unnaturally perceived to be almost a branch of science. Now, in its relative 
maturity, it is, or at least it should be, relegated to what it really is: a very 
valuable technique for drawing scientific inference. The basic physics of seismic 
oscillations is not new, is fundamentally quite simple, and is well understood; 
and the principles of inference from the observations should be straightforward 
to comprehend, even though the technicalities of putting them into practice may 
for some seem to be rather complicated. Therefore the robust raw conclusions 
are genuinely secure, more so than the broader issues to which they are intended 
to be applied. However, it is incumbent upon us to make the distinction between 
the inferences that really have been drawn reliably and the further wider infer- 
ences that might subsequently be, or have been, drawn, often with the aid of 
supplementary, possibly less secure, maybe non-seismic, information, and even, 
maybe, (sometimes unstated) surmise. Only if such distinction is made clear can 
the contributions of our subject to science be reaped to the full. Unfortunately, 
that aspiration has not always been achieved in the past, and misinformation 
has sometimes sullied the waters. 

It takes only a brief scrutiny of the equations describing the structure and 
dynamical evolution of the Sun (it is not quite so brief to derive them) and the 
equations governing the low-amplitude seismic modes of oscillation to appreciate 
what broadly can, at least in principle, be reliably inferred. Anything further 
must depend on other criteria, such as general physical argument beyond seis- 
mology, traditional astronomical observation, or even prejudice. It is obligatory 
to be explicit about how such additional constraints are applied. The subject has 
advanced to a new level of sophistication; we are now trying to probe seismically 
(and otherwise) almost inaccessible aspects of the physics of the Sun, and the 
techniques for unravelling them are becoming more and more intricate, beyond 
the point at which most scientists wish to tread. There must necessarily be an 
increased trust in our findings, and it is our responsibility not to betray it. Many 
of the broader scientific community want to use our results in their research; for 
that they need to know not only the limitations of our inferences, and the caveats 
upon which they are based, but also which aspects of what we seismologists tell 
them can really be trusted. 

Much of the emphasis of Solar Dynamics Observatory seismology concerns the 
workings of the convection zone. We want to know what controls the solar cycle, 
how magnetic field is amplified, modulated, and then suppressed, how sunspots 
are formed and destroyed - and what determines their lifespan. We want to 
know the geometry of at least the larger scales of convective motion, and how, 
beneath the seen superficial layers of the Sun, the processes that control the 
total radiative output are modulated. At least some of us want to understand 
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how all these matters influence our procedures for inferring the gross properties 
of the Sun, and how they impinge on our broader ideas of the evolution of 
the Sun in particular, and of stars in general. Addressing such delicate issues 
with confldence may now seem an almost impossible task to us older scientists 
who have lived through the years of stumbling through the darkness, having 
finally emerged to bathe in the secure light illuminating the minute arena of 
knowledge that we have been instrumental in uncovering. It is now up to the 
younger community to proceed likewise: to grasp at the edge of our perception 
with initially insecure ideas, fully appreciating the uncertainty, of course; then 
moulding and strengthening them into a new body of secure scientific knowledge. 



2. Introduction 



Once the p otential of solar oscillations to map the interior of the Sun was 
recognized ( Christensen-Dalsgaard and Gougb . 19761 ) there were two obvious 



serious issues that immediately appeared accessible to resolution: the spherically 
symmetric component of the hydrostatic stratification and the internal angular 
velocity. The first of these was needed for investigating what has been called the 
solar neutrino problem; the second concerned the centrifugally induced oblate- 
ness of the Sun's gravitational equipotentials, and how that impinged upon an 
important test of theories of gravity. General Relativity in particular, via the 
precession of planetary orbits. These were the two most widely discussed issues 
in heliophysics at the time. 

The st r atifica tion was originally addressed with the help of theoretical models. 
Deubner d 19751) had published the first well resolved k-uj spectrum, and Ando 



and Osaki ((1975I ) had shown that theoretical eigenfreguencies of p modes trapped 
in the outer layers of a solar model envelope were in quite good, but not perfect, 
agreement with Deubner's observations. What was required for bringing the 
theory more closely into line with observation was first estimated from the prop- 
erties of the eigenfrequencics of a simple polytropic representation of t he outer 



layer s of the convection zone as a means of calibrating solar models (jGoueh 



19771 ): the conclusion was that the convection zone must be about 200 Mm 
deep, some 50 Mm deeper than the favoured value of the time. Basically, that 
conclusion was drawn from relating theoretical eigenfrequencies to the jump 
in the adiabatic "constant" p/p'^^, which is closely related to specific entropy, 
across the thin superadiabatic convective boundary layer. It was subsequently 
suppo rted by more realistic, numerical, computations by lUlrich and Rhodes! 



(jl977f ). The principal implication of that result was that, according to complete 
solar models, a deeper convection zone implied greater helium and heavy-element 
abundances, a hotter, more centrally condensed, core, and a higher neutriiio flux . 



thereby exacerbating the solar neutrino problem {e.g. Abraham and Ibenl . 1971 



Bahcall and Ulrichl .[ I971I ). Indeed, one is tempted to speculate that the modellers 



in the past had adjusted the defining parameters of their models to minimize 
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the theoretical neutrino flux, and that that had prejudiced Ando and Osaki's 
calculations, although without repeating the calculations oneself (and maybe 
even if one did) one cannot be sure. 

There was a great deal of healthy mistrust in the models at that time. The 
concern was that a seismological measurement - I should really say estimate - 
of merely the upper boundary layer of the convection zone, extending only a 
minute fraction of the solar radius beneath the photosphere, could hardly be a 
robust indicator of conditions inside the energy-generating c ore. Therefore, as 
soon as lo w-degree data became available (IClaverie et al. ^ 1980 ^Grec, Fossat, and 
Pomerantz, 19801 : Fossat. Grec. and Pomerantzl . 198lf ) it became possible to use 
seismic indica tors of mo re global propert ies; first the so - called large frequency 
separation ((Christensen-Dalsgaard, Gough. and Morganl. ll979aHbl: Christensen- 
Dalsgaard and Gough, 198ot ). which measur es the sound travel time from the 
centre of the Sun to the seismic surfac43 {e.g. 



Vandakuroyl. 11967 



Tassoull . [l980l : 



iGoughl . Il986l ). and later the small frequency sepa ration (Gough . 1983c ). which 
is a direct indicator of conditions in the core {e.g. Gough . 1983a . 19861 ). Unfor- 
tunately it was not possible to make all of the theoretical frequencies agree with 
the data, as was clear from an indiscrimina nt attempt to fit only the low-degree 
frequ encies with whole-disc observations (jChristensen-Dalsgaard and Goughl . 
19811 ): two solar models seemed to be favoured, one with a low initial helium 
abundance: Yq = 0.18 (coupled with a correspondingly low heavy-element abun- 
dance Z = 0.003 and a low neutrino flux, although not low enough to reproduce 
the neutrino detection rate), the other large: Yq = 0.27 and Zq — 0.026. The 
high-1" model fitted rather better, and one was tempted to prefer it, especially 
because its high-degree mode frequencies were closer to Deubner's observations. 
Moreover, the low-F model had a helium abundance below what was thought to 
have been produced in the Big Bang (c/. footnote 8), which would call for some 
contrived explaining. 



consider the seismic surface r = R of the Sun (assumed here to be spherically symmetrical) 
to be the radius at which c^, regarded as a function of r, or c, regarded as a function of 
acoustic radius T(r) = / c~^dr - both of which are cl ose to being linear functions in the outer 
adiab aticallv stratified layers of the convection zone llBalmforth and Go ugh. 1 990t Lopes and 
Gough^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^Fchristensen-Dalsgaard 
it lies about 1000 km above the photosphere, the precise value depending on 
exactly how the extrapolation is carried out. There is nothing special about the structure of 
the actual atmosphere in its vicinity, which is well inside the outer evanescent zone of most 
of the seismic modes and therefore has little significant influence on the dynamics. Instead, 
it acts simply as a (virtual) singularity in the acoustic wave equation, providing a convenient 
parametrization of conditions (well below the photosphere) in the vicinity of the upper turning 
points of the modes. Put another way, it provides a convenient fiducial location with respect to 
which the acoustic phase in the propagating zone beneath is related. Unlike the photosphere, 
which has no acoustic significance, it shares a relation with the deeper solar interior that 
is robust, and is insensitive to the non-seismic, thermal and radiative, properties of the outer 
convective boun dary layer , whose stru c ture changes with the solar cycle llAntia and BasuLl2004l : 
iDziembo wsici and Goo del l2004l 120051; iLefebvre and Kosovichevl. I2005t Lefebvre, Kosovichev , 
and Rozclot, 20071. In contrast to other, non-seismic, radii (cf. iBahcall and Ulrichl. Il988l) . 
it provides a stable outer limit to the effective total acoustic-radius integral t{R), which 
determines the large frequency separation; in Model S it is some 200 seconds or so greater 
than the actual acoustic radius of the photosphere. 
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The two models required different identifications of tlie orders [n] of tlie 
modes; unfortunately the orders were so high, and the fit s o poor, that reliable 



extrap olation to n = 1 was not possible. It was not until iDuvall and Harvey 



(jl983f ) observed the frequencies of modes of intermediate degree that a secure 



connec tion between modes of low and high degree could be made (jGoueh 



1983B). The orders of the latter are determinable because the frequencies of 



high-degree f modes, with which is associa ted n = 0, are essentially independent 



of the structure of the Sun (jGoughl . Il982al ) . The high-F alternative was thereby 
confirmed. 

Thus we had learned that the solar neutrino problem was almost certainly 
not resolvable by adjusting solar models, and must be a matter for nuclear or 
particle physics. 

Had we really learned that? Perhaps not yet. At the time there were still very 
serious doubts about the solar models, for they depended on many unproven as- 
sumptions, some of which are listed in Table 1. There had already been, and there 
were yet to be, many models computed in which some of these assumptions were 
relaxed in the hope of yielding lower neutrino fluxes (although none provided a 
satisfactory reconciliation of theory with observation). At the very least, a secure 
representation of the stratification throughout the Sun was surely required. That 
was soon to be provided from inversion a.nalyse s of frequency data from Duvall 
and Harvey ( Christeiisen-Dalsgaard et all. 19851: Christensen-Dalsgaard, G ough, 



and Thompson, 19891 : Dziembowski. Pamvatnvkh. and Sienkiewicz . 1990l ). fro 



3m 



which it was possible to infer the (spherically averaged) sound speed [c(r)] 
throughout almost all of the Sun. Details of the core were not yet within reach 
(and even today there is considerable uncertainty), but elsewhere the sound 
speed was essentially the same as that of Christensen-Dalsgaard and Gough's 
(1981) high-Y model. It was also possible to sec the base of the convection 
zone. I recall one Friday morning (the second Friday of January, 1984, a day on 
which I was due to deliver a lecture on my findings to the Royal Astronomical 
Society); it was about 4:30 in the morning - one could compute seriously only 
at night in those days - when I obtained my first plot of c^(r) that extended 
beneath the convection zone; I plotted and not c because it is related closely 
to temperature T - the equation of state in the solar interior is reasonably well 
approximated by the perfect-gas law, for which oc T/ fx, where is the "mean 
molecular mass" . A more modern (and rather better) version of what I obtained 
is illustrated in Figure [J. The location of the base of the convection zone is 



^Included for comparison is one o f the models with low Y that was used in the original cali- 
bration with low-degree modes bv lChristensen-Dalsgaard and Goughl | |1981D . The qualitative 
differences between the models can be appreciated by realizing first that the radiative envelopes 
are roughly polytropic (with index 3.5), and that it is adequate to approximate the equation 
of state by the perfect-gas law. Then it is clear that the magnitude of T{r) (and p) must be 
greater in the higher-Y model, because the total energy generation rate, which i s an increasing 
function of X, p, and T, is the same for the two models. Polytropic scaling (e. a. iGoughl . Il990f) 
indicates that Z is a steeply decreasing function of X, so Z is much greater in the higher- Y 
model. In addition, the polytropic radius scale is greater for the higher-V model, as is evident 
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Figure 1. Upper panel: The continuous curve extending over the entire range of r/i? is 
the S(juareo^^hc_soun^_s£ee^Jc^j_j^^he_sto 

et al. 1119961) . which was computed with initial hydrogen and heavy-element abundances 
X(j = 0.7091, Zq = 0.0196. The square of the sound speed in the Sun is also plotted (where 
r/i? > 0.1) as a continuous curve; the difference between the two can be barely discerned by the 
variation in the thickness of the apparently single curve (but see Figures |3] and |4]l . The dip in 
at the centre of the Sun is a result of the augmentation of by nuclear transmutation, which 
increases with time on the main sequence (and which can therefore be used as a diagnostic of 
main-sequence age). The dashed curve is on the zero-age main sequence, and, except near the 
surface where abundant elements undergo ionization, is (locally) proportional to temperature. 
The dotted curve is in a model with an initial heavy-element abundance Z = 0.001, 
but continuously contaminated at the surface at such a constant rate as to have a current 
heavy-element abundance_^^_^_0j02jr^^lie_^onvectMi^_z^ 

Gough, and Morgan TTgTgalK LowerTaneir^^mDcrat^ two present-day solar 

models. 



evident as a near discontinuity of the second derivative, as one can see more 



from the Figure ^ by imagining an extrapolation of the functional form of T(r) outwards 
from the radiative zone. Consequently the convection zone, which steepens the gradient, has 
more truncating from the radiative structure to perform in order to maintain the observed 
photospheric radius, and is therefore deeper. An additional scaling in magnitude is required to 
convert T to oc T / fi, raising the dotted curve (corresponding to the model with the lower 
Y) relative to the continuous curve. It also depresses both curves near the centre of the star, in 
the energy-generating core where fi has been augmented by nuclear transmutation, providing 
a diagnostic of main-sequence age. 
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easily by holding the page almost in one's line of sight and looking along the 
curve. I confirmed the location by having a research student, the only person 
left in the computer room at that hour, repeat the exercise. I went home, slept 
a couple o f hours, and th en took a train to London. When I showed the plot at 
the RAS (ICoughl Il984ah and pointed out the discontinuity, the audience was 



incredulous, even though I stood almost in the plane of the screen to reassure 
them. I then recounted the confirmation by the research student, and there was a 
sudden release of tension in the audience when I told them that the student was 
a cosmologist, for then they appreciated that such a student must undoubtedly 
have provided an unbiased opinion on such a matter. 

T he result was compared with theoretical solar models f Christenscn-Dalsgaardj 
et al.. \l98^ . and revealed a characteristic discrepancy (much larger than that 



evident in the more modern comparisons, such as those illustrated in Figures [3] 
and H]) immediately beneath the base of the convection zone. It was suggested 
that the discrepancy could have been caused by an error in the opacity (values 
of opacity in those days were available in tables provided by the Los Alamos 
National Laboratory in the USA). Such sugg estions, i n a r ather broad sense, 
had been made in the past, for example by ISimon (Il982l) in response to a 



failure to reproduce simultaneously observations relating to stellar evolution 
and observations of pulsational characteristics of classical variable stars. But 
here one was able to state quite precisely the thermodynamical conditions under 
which the opacity was in error, and the approximate magnitude (and sign) of the 
error. After some persuasion, Carlos Iglesias and Forrest Rogers at Livermore 
computed by moonlight the opacity for a few judiciously chosen values of the 
state variables p and T (and a plausibly appropriate chemical composition) using 
an independently written computer programme; they not only confirmed the 
helioseismic inference but also showed that the discrepancy was even greater 
at lower temperatures - of no concern for the Sun because it is well inside 
the adiabatically stratified region of the convection zone where radiative energy 
transfer is negligible, but of great significance in reconciling with observation the 
theory of some cl asses of int rinsicallv variable stars such as Cephei a nd slowly 
pulsating B stars dCox et ali[l 992: Ki riakidis. El Bid, and Glatzeill992l: Moskalik 



and Dziembowski. Il992l : Dziembowski and" Pamvatnvkhl. 19931: Dziembowski. 



Moskalik, and Pamv atnvkh. Il993r 19941) and double-mode Cepheids f Moskalik. 



Buchler, and Marom, 1992f ). Largelv as a result of efforts bv Werner Dappen, 



scientists who were involved with the Los Alamos opacity computations were 
brought into the same room as the Livermore scientists, and errors in the Los 
Alamos calculations were identified. This led to the Lawrence Livermore National 
Laboratory provi ding the funds for m ore extensive opaci ty computations for 
astrqphysical use (llglesias. Rogers, and Wilson . ,1990. : Jglesias and Roger l ll99l[ 



■ l2002f ). Thus we see the first example of helioseismoL 



ogy contributing directly to microscopic physics (if one considers the neutrino 
issue to be indirect). There have been further contributions, but I postpone 
discussion of those until later. 

I emphasize that an important consequence of the sound-speed inversion was 
that it convincingly ruled out the \ow-Y models, reaffirming that the cause of 
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the neutrino deficit was not in solar modelling. That conclusion was not fully 
appreciated by many of the solar modellers of the time, for they had not yet 
understood the power of helioseismological analysis. 

I turn now to the angular velocity [17], which is measured by the odd (with 
respect to azimutlial order m) component of degeneracy splitting caused by 
advection and Coriolis acceleration. Strictly speaking, one needs the sound-speed 
inversion first, to establish the hydrostatic stratification with respect to which 
the splitting kernels are computed. However, the angular velocity was actually 
first obtained by using a theo retical solar model before the full sound-speed 
inversion had been carried out ( Duvall et al. . 19841 ). partly, perhaps, because it 



was easier to perform with the data in hand, yet also with the knowledge that 
the error from which that procedure suffered was much less than that resulting 
from the errors in the measurements of the degeneracy splitting. The result was 
a big surprise: 



There had been many prior discussions of how much faster than the pho- 
tosphere the solar core must be rotating (due to spin-down), the most widely 
publicized being the discussions resultin g from the measurements of t he a' 
parent oblateness of the solar surface by Dicke and Goldenberg ( 1967 , Il97 
The oblateness measu rements had been made in the hope of confirming Dicke's 
theory of gravitation ( Brans and Dickj . ll961 lR: put naively, in the Brans-Dicke 



^The gravitational attraction associated with the energy density ~GM/r of the gravitational 
field surrounding the Sun, absent in Newton's theory, causes the total gravitational attraction 
to increase: very roughly speaking, as a result of energy conservation the apparent gravitational 
mass of a planet at distance r from the Sun, in a fiat representation of space, is augmented 
by approximately GM/c^r per unit mass above what it would have appeared to have been 
at infinity; M is the mass of the Sun, and here c is the speed of light. Similarly, the energy, 
hence the frequency, of a photon is multiplied by a factor F = 1 -|- GM/c^t - that causes 
the familiar gravitational redshift. Consequently, the orbit equation is modified simply by 
multiplying the Newtonian gravitational force on the planet by F^. After linearization and 
rewriting M in terms of the orbital specific angular momentum h = ^J\GMr), valid for nearly 
circular orbits, the effective attractive force becomes — (1 -|- 3/i^/c^r^)GM/r^; it increases with 
increasing proximity more rapidly than Newton's inverse square. It is easy to see also that the 
gravitational field in the equatorial plane of a rotating (oblate) axisymmetric self-gravitating 
body (like the Sun) also increases with decreasing distance faster than the field around a 
corresponding spherically symmetrical body: the act of flattening a spherical body takes equal 
amounts of material from the poles towards the near and the far sides of the equator, the 
increase in gravitational attraction by closer nearside matter exceeding the lesser decrease by 
farside matter. Therefore the net gravitational attraction is increased, by an amount which 
inceases as r decreases and the shape of the Sun becomes more apparent. The force is given 
approximately by —(1 -|- ^ J2iJ^/r^)GM/r-^; J2 is the quadrupole moment. In both cases, 
therefore, a planet is drawn towards the Sun increasingly strongly with increasing proximity 
than it would have been in an inverse-square field. Conserving its angular momentum, it is 
thereby caused to rotate through a greater angle near perihelion because its orbital angular 
velocity is augmented, distorting an otherwise approximately elliptical bound Newtonian orbit 
in such a manner as to appear to make it simply precess in the same direction as the angular 
velocity of the planet. (Near aphelion the oppositely directed contribution to the precession 
is lesser, therefore too small to annul the contribution from near perihelion.) Use of planetary 
(or spacecraft) orbital precession rates to calibrate the relativistically induced deviation from 
the inverse-square gravitational field surrounding the Sun therefore requires one to know the 
contribution from the distortion from spherical symmetry of the mass distribution in the Sun, 
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theory Newton's gravitational constant [G] was regarded as a field satisfying 
a wave equation that couples to the matter, thereby relaxing the gravitational 
force field from the relative rigidity that is imposed when G is held constant, and 
hence reducing the rate of precession of planetary orbits from the value predicted 
by General Relativity; the fact that General Relativity predicts essentially the 
correct precession with a spherical Sun was regarded by Dicke as a fortuitous 
coincidence, and that in reality the shortfall predicted by his theory (which 
is not an absolute prediction, but depends on an unknown coupling constant) 
is made up by the oblateness of the Sun's gravitational field produced by a 
rapidly rotating interior. The oblateness measurements would calibrate the cou- 
pling constant. A debate ensued concerning two matters: the relation between 
the oblateness in surface brigh tness and the oblateness of the gravitational 
equipotentials {e.g. Dicke . 1970l ). and the fluid dynamics of the solar interior 
conc erning angular-momentum tra.nsfer f r om centre to surface during spin-down 



Howard. Moore, and Spiegell . 119671 : iDickd . 119671 : iBretherton and Spiegel 



1968|). I shall enlarge on neither here, desp ite the intrinsic int erest of each, 
because even the early seismological analysis (jPuvall et all Il984l ) revealed that 
the Sun is unexpectedly rotating almost uniformly throughout its interior, except 
possibly in what might be its almost inaccessible core (Figure ^ . Subsequent 
observations of low-degree modes revealed that the core rotates significantly n o 
more rapidly than the envelope ( Elsworth et al . 19951 : IChaplin et al . 1999^ . 
From knowledge of fl and the density and pressure stratification, the multipole 
moments J2k of the gravita tional equipotentials can easily be computed {e.g. 
Goughl . 1981 : Piipers . 1998[ ). I maintain that this is the most accurate way to 
determine J2k, notwithstanding some claims in the literature to the contrary. 
The reason is partly that other methods involve relating surface brightness 
asphericity with gravity, which is not well understood, and partly because those 
measures rely on only the small quadratic centrifugal force rather than the more 
robust, and unambiguous, seismic effects of advection and, to a much lesser 
degree, Coriolis acceleration, both of which are predominantly linear in 17 and 
are therefore much larger than the centrifugal term; and, moreover, they do not 
depend on relating brightness to gravity. Moreover, the rotational distortion of 
the isobaric surfaces in the vicinity of the photosphere is dominated by the direct 
effect of the local centrifugal force, and exceeds the gravitational distortion by 
a factor of more than 20; determining the latter from the apparent shape of the 
solar disc therefore requires a subtraction of two measurements that differ by less 
than a mere 5 per cent. Even the mo st precise limb-shape measurements { e.g. 
Fivian et aR l2008l : iKuhn et all |2012[ ) are far from that goal (jCoughl . 12012^ ). 



The outcome is consistent with General Relativity; indeed, the minisc ule con- 
tribution that J2 does make to the precession of the orbit of Mercury (e.o.lGough , 
1982bl : iDuvall ei all . 1 19841 : iBrown et all . 1 19891 : lAntia. Chitre. and Goughl . I2OO8I ) 
brought General Relativity into somewhat closer agreement with observation 
than assuming the Sun to be spherically symmetrical. The influence of higher 



such as is produced by centrifugal accelerati on due to rotat ion. The precession rate is most 
easily calculated by perturbation theory (e.g. iRamsevl [T937D . 
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Figure 2. Optimally localized averages of the North -South symmetric component of the 
angular velocity [Q] of the Sun, depicted as a contour plot in a quadrant. Some contours arc 
labelled [nHz], and, for clarity, every fifth contour is drawn bold; the contour separation is 
lOnHz. The outer quarter circle denotes the surface of the Sun, the dashed quarter circle 
indicates the base of convection zone, and the tick marks at the Sun's surface are drawn at 
latitudes 15, 30, 45, 60, and 75°. The Equator is the horizontal (relative to the page) axis 
and the Pole the vertical, each labelled with values of x = r / R. The shaded area indicates 
the region in the S un where no reliable inference could be made from the data available (from 
ISchou et a71ll998D . There has been little significant i mprovement of this in f erence, save for 
a dem onstration that the core is not rotating rapidly llElsworth et al\ . 119951 ; IChaplin et all 
Il999t): however, there have been studies of temporal variation associate d with the solar cycle 
{e.g. IVorontsov et all \200i iBasu and Antial 120031; iHowe et all l2006allbt Antia. Chitre, and 
Gough. \200K 

moments is utterly negligible. The improvement is marginal, but with future, 
more delicate, orbital analyses the helioseismological determination of J2 will 
evidently take on a crucial role. The almost uniform rotation of the radiative 
interior of the Sun begs the question of how that can be, particularly because 
it interfaces with a differentially rotating convection zone. I shall return to that 
matter later. 



3. Understanding Seismic Variables 

It is extremely important to understand what is actually inferred from seis- 
mology, for only then can one draw reliable conclusions concerning the Sun. 
Seismic modes are essentially adiabatic. They result from forces - pressure and 
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gravity predominantly - acting on matter with inertia. The dynamics there- 
fore concerns only pressure [p] and matter density [p] which are related via 
gravity through hydrostatic equilibrium, and the relation between perturbations 
to them under adiabatic change, which is characterized by the first adiabatic 
exponent [71 = (91np/(?lnp)^, the partial thermodynamic derivative being 
taken at constant specific entropy [s] ]. The restoring force of acoustic modes 
(p modes) is principally pressure, that of gravity modes (g modes, including 
the fundamental f modes) is buoyancy. So the modes can provide information 
directly about only p,p, and 71 (and, of course, any function of them). These 
are the basic seismic variables. I should acknowledge that any magnetic field [B] 
that is present also contributes to the dynamics, and therefore is also a seismic 
variable. Unfortunately, its effect on the frequencies of the modes appears to 
be indistinguishable from that of an appropriate variation in sound speed [c] - 
also a sei smic variable because = lipl P - which makes it difficult to unravel 
the two ( Zweibel and Goughi . [Hii). In principle one might be able to do so 



from the eigenfunctions, which are configured differently by c and B, although 
supplementation by non-seismic arguments is likely to be more productive. 

Leaving B aside for the time being, it should be appreciated that the relation 
between p and p through hydrostatics does not depend directly on 71, so p 
or p (or any function of only them) and 71 are structurally independent: it 
follows that although the relation between them, which physically is given by the 
equation of state, can be determined by seismology along the thermodynamic p- 
p path through the Sun, there is no remaining redundancy, so the veracity of the 
equation of state itself cannot therefore be probed by seismology alone; in order 
to investigate the equation of state, supplementary, non-seismic, information is 
required. 

There are two approaches that one can take for drawing seismic inference. One 
is to adopt a parametrized model of the Sun, or some aspect of it, and from it cal- 
culate whatever seismic properties one wishes to compare with observation. The 
comparison calibrates the controlling parameters. That was the procedure that I 
described in my introduction for first estimating the depth of the convection zone, 
and consequently the helium abundance. The other approach is to ignore explicit 
models (almost) entirely, and combine the data in such a manner as to isolate 
certain properties of the seismic structure. Typically that involves adopting, at 
least at first, a reference model of the Sun, which it is hoped is sufficiently close 
to the Sun for linearization in the small differences from it to be more-or-less 
valid; that simplifies the analysis enormously. Subsequent iteration can usually 
remove the dependence on the reference model. Nonlinear asymptotic methods 
have also been used successfully to yield approximate inferences without recourse 
to a reference model at all. The most common procedures that have been used 
to date are aimed at obtaining easily interpretable representations of the basic 
seismic variables. They are commonly called inversions. 

When linearization about a reference model is carried out, the frequency 
differences between the Sun and the model can each be expressed as a sum 
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of spatial averages of independent seismic variables. A prudent procedure is 
then to seek suitable linear combinations of those averages which more easily 
inform one of whatever question one has chosen to pose. If it is the value of a 
particular seismic variable that one wants to investigate, then the most easily 
interpretable data combinations are those that represent averages of that variable 
with weighting functions (called kernels in this context) that are highly localized 
with hardly any sidelobes, and which at the same time are (almost) independent 
of any other variable, for then the result can be thought of as a blurred view of the 
variable of interest. Consequently one tries to tailor maximally localized kernels, 
although the attempted maximization must be moderated by a requirement 
that the interference from extraneous seismic variables be kept low and that the 
influence of data errors not be excessive. On the whole, increasing localization 
requires data combinations with coefficients of greater and greater magnitude, 
which increases the influence of random data errors. The procedure to construct 
those kernels is now called optimally localized averaging (OLA). Just how low 
one demands the interference from other variables to be, and how much the 
influence of data errors one considers to be acceptable, is a matter of personal 
choice, which probably explains why the errors in the averages displayed in 
Figures [3] and U] differ. This ki nd of tailoring was the first to be advocated for 
helioseismology (jGoughl . 119781 ). in contrast to the almost universal opinion of 
geoseismologists of the time who objected to the procedure on the ground that 
one might be tempted simply to draw a curve through the averages and mistake 
that curve for the actual variable, rather than the average that it actually is. 
We helioseismolgists usually do not make that mistake. We explicitly draw a 
sequence of crosses, as in Figure [3l the horizontal components representing the 
averaging widths, the vertical bars representing the propagated standard errors 
in the data. Some geoseismologists are now changing their tune, and are following 
suit. 



An alternative procedure, common in gcoseismology, is to seek a putative 
"solution" that reproduces the data. Typically, one first expresses the seismic 
variables as linear combinations of a prescribed set of basis functions, and then 
one chooses the coefficients in each combination to match the data the most 
closely, again moderated (regularized) to avoid excessive error magnification. 
The reason geoseismologists used to prefer that approach is that it might yield 
a curve that actually fits the data within the estimated errors, which OLA 
might not (and is not explicitly designed to). If it is a linearized perturbation 
to a reference model that is being fitted to the data, the outcome can also 
be expressed as averages with localized components. However, the averaging 
kernels often have severe sidelobes too, usually near the surface of the Sun, 
which renders interpretation more difficult. Moreover, the "solution" , if it exists, 
is just one of infinitely many that satisfy the constraints imposed by the data. 
Usually, error-correlation information about the raw data is unavailable, and the 
data are usually fit by (regularized) least squares (RLS), regularization being 
accomplished typically by imposing some criterion of smoothness on the ground 
that the data provide only finite spatial resolution. 
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Figure 3. Optimally localized averages of the relative differences of the squared sound speed 
[c^] and the density [p] in the Su n from those in Mo d el S of Christensen-Dalsgaard et al. 
(1996), computed by M. Takata jTakata and Goughl, 1200 Ih from MDI 360-day data and 
plotted against the centres [x = f/R] of the averaging kernels [A(x;a;)] (which here resemble 
Gaussian functions), defined hy x = J xA^dx / J A^dx. The length of each horizontal bar is 
twice the spread s of the corresponding averaging kernel, defined as s = 12 J(a; — x)^A'^dx 
- an averaging kernel A that is well represented by a Gaussian function of variance has 
spread approximately 1.7A ^ 0.72 FWHM; were it to be a top-hat function, its spread would 
be the full width, which is why s has been so defined. The vertical bars extend to ±1 standard 
deviation of the inversion errors, computed from the frequency errors quoted by the observers 
assuming them to^^stati^ticaU^mdegendent^tl^^ 
Sekii, and Stark. 1199631 1 Howe and ThomDsonl.ll996D . 



If two seismic variables are involved, as is the case for structure inversions (in 
the absence of rotation and a magnetic field), then two separate complementary 
inversions must be carried out for OLA, each requiring the suppression of the 
contribution to the data combination from the variable not being sought. With 
the help of some regularizing assumption, estimates of the seismic variables can 
then be made from their averages to enable one to gauge, for each seismic average, 
the contribution to the corresponding data combination from the other seismic 
variable; in most of the published "inversions" , that contribution is ignored. In 
the case of RLS frequency fitting, the two seismic variables can be represented 
simultaneously. 

I do not here go into how the inversions arc carried out. That is not necessary 
for understanding the results, provided that adequate information about the av- 
eraging kernels is given. Unfortunately, that information is not always provided. 
Precisely how the centres of the optimally localized kernels were determined is 
rarely stated, although if the kernels are very narrow it doesn't much matter, un- 
less there are large distant sidelobes. Sometimes authors use merely the locations 
where they tried to centre the kernels, rather than where they actually succeeded 
in centring them, which requires conjecture on the part of the reader relating to 
the information contained in the data set employed and the proficiency of the 
author in extracting it. One useful rule of thumb for helping to guess what might 
have been plotted is to recognize that with currently available frequency data 
well localized kernels without substantial sidelobes cannot be centred closer to 
the centre of the Sun than t/Rk, 0.05. 
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Figure 4. Several sets of optimally localized averages of the relative differenee between the 
squared sound speed in the Sun and that in Model S of Christensen-Dalsgaard et al. (1996) 
computed from the same frequency data set, with kernels obtained by different inverters and 
by ignoring the contamination from other seismic variables. The vertical bars extend to ± one 
standard deviation of the errors computed from the frequency errors quoted by the observers. 
[The abrupt deviat ion close to the surface is a characteristic of having misrepresented the 
acoustic radius fcf. iTakata and Goughl . I2003D .] 

Inversion is not a well prescribed procedure. It is an art. And it can lead to 
different representations of the information, as is illustrated in Figure HI which 
depicts different optimally localized sound-speed averages all derived from the 
same data. It should therefore be recognized that the differences between the 
averages plotted in the figure are not necessarily an indication of inversion error; 
instead, they can result from the differences in the averaging kernels selected. 
It is evident, therefore, that merely offering some measure of the widths of the 
kernels, which is absent from Figured is insufficient for appreciating the results 
fully. 

The objectives of most of our investigations do not concern seismic variables 
alone. I have already mentioned issues concerning the solar neutrino flux, and 
also the helium abundance [Y] and the relation between its initial value [Yq] 
and Big-Bang nucleosynthesis. Investigation of these non-seismic quantities using 
seismology necessarily involves relating them to seismic variables via a theoretical 
model, which itself depends on the assumptions upon which that model has been 
built. That statement may seem obvious. But those not in our subject do not 
always appreciate just which of those assumptions are important. It therefore 
does no harm to state them. Evidently, which of the assumptions are the most 
important depends on the matter in hand. 

I conclude this part of my discussion with another simple point, which is also 
not obvious to everybody: the dominant physics of helioseismology is extremely 
simple; it is simply the physics of the propagation and interference of well un- 
derstood waves. Therefore, provided one is scrupulous in presenting the results 
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with due care and attention, the direct conclusions cannot be questioned. When 
an observation appears to be in conflict with our seismological knowledge, as 
superficially seems to be the case with modern spectroscopic determinations of 
the photospheric heavy-element abundances, for example - an issue to which I 
shall turn my attention below - con trary to the opinion of some commentators 



{e.g. iGuzik, Watson, and Coxj . 120061 ) it is not the seismology itself that is to be 
challengeco. I hasten to add that there certainly are issues within seismology 
that have not been resolved, usually because a suitable way to analyse the data 
has not yet been found (possibly because the data do not even contain that 
information in a form that permits it to be readily extracted, possibly because 
it is impossible to do so, possibly because the extraction procedure itself still 
seems to be beyond our capabilities), but that is a different matter. 

One further matter concerns the manner in which wc use a diagnostic. That 
typically depends on the issue one wishes to address, and its relation to that 
issue may itself be subject to some doubt, perhaps due to untested assumptions 
in modelling. This means that the accuracy of the result may be much less than 
the precision, perhaps very much less. It is important to recognize the difference 
between the two. Precision can usually be estimated well from the precision of 
the data - although in order to trust the outcome one has to trust the estimated 
errors in the data. Accuracy requires recognizing and assessing the influence 
of the assumptions - inaccurate assumptions can lead to a dispersion amongst 
the outcomes in any in vestigat i on that prudently utilizes a variety of different 
analysis techniques (c/. iGoueh . 2012al ). Therefore precision is never judged to 



be lesser when fewer analysis procedures are considered. Unfortunately, greater 
precision is often mistaken for greater accuracy. 



4. Standard Solar Models 



Standard solar models arc constructed usually with the most sophisticated mi- 
crophysics available, yet with the most primitive macrophysics - fluid dynamical 
processes are ignored wherever possible, probably because they are too difficult 
to model in an agreed standardized manner. Rotation and magnetic fields are 
normally ignored too, so the star is spherically symmetrical. The star is evolved 
hydrostatically either from somewhere on the Hayashi track, gravitationally 
contracting until it reaches the main sequence, whence radiative energy loss 
is balanced almost exactly by nuclear energy generated in the core; or it is 
evolved from an estimate of the zero-age main-sequence structure, which is 
obtained by equating radiant energy loss with nuclear generation in a completely 



*I am assuming that the seismic data proccssers have judged their errors correctly, and that 
the error correlations in the frequency data sets, normally ignored, are not unduly severe. I 
am also assuming that, where it is appropriate, asphericity of equilibrium structure is taken 
correctly into account. 
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Figure 5. The dashed Hnc is the relative difference between the sound speed i n a standard 
solar model computed by Bahcall and Pinsonncault and that infer red by Basu llBasu et al\ . 
[19970 froni a combinat ion of osc illation fre quency data obtained bv lChaplin et alTi^gg^ and 
iTom czvk. Schou. and Thomgson|j^995a); Tgmcz^l^_e|_aiJJl995b|]_(fron^ B ahcall , Pinsonneault, 
and Basu, ,2001}. The arrow labelled ^Bc etc. is simply a very rough estimate of the sound-speed 
discrepancy inatheoretica^_sola^niode^_adji^^ Basu, and 

Pinsonneault ||1998D to yield a neutrino flux concomitant with that measured by the gallium 
detector at Gran Sasso l lHampel et a;.LI1996D . 



homogeneous star. The second approach is not quite correct, because there was 
some nuclear transmutation prior to arrival on the main sequence, most notably 
conversion of '^He to *He, which temporarily slows down the contraction (but is 
InsufRcient to halt it); however, that phase has little impact on the subsequent 
evolution of the star. Although the transition from gravitational to nuclear 
energy release is smooth, during subsequent evolution on the main sequence 



the central hydrogen abundance [Xc] declines almost linearly with time (jGouehl . 



19951 ) , so backwards extrapolation of to the initial hydrogen abundance [Xq] 



provides a useful fiducial origin of main-sequence age. 
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Tabic 1: Standard model assumptions 



Initially uniform chemical composition 
Spherical symmetry 
Hydrostatic equilibrium 

Simple description of energy transport, using mixing-length theory in convection 
zones, and a simplified treatment of radiative transfer in the atmosphere 

No serious internal mixing of chemical species except in convection zones, 
so no tachocline, yet gravitational settling and chemical diffusion 

Rotation dynamically negligible 

Maxwell stresses negligible {i.e. no magnetism) 

No mass loss 

No accretion 

No large-scale instability 

No (nonlinear) transport by waves 

Thermal balance (almost) 

Equation of state known 

Nuclear reaction rates (reaction cross-sections) known 
Opacity formula known 



Evolution on the main sequence is hydrostatic, and essentially in thermal 
balance, energy being transported from centre to surface by optically dense 
radiative transfer wherever the stratification is convectively stable: I present the 
transfer equation in the next section. Nuclei undergoing thermonuclear reactions 
are usually considered to be screened by e lectrons accordi ng to the classical, non- 
relativistic, Debye-Hiickel treatment by ISalpeteij (Il954l). although the validity 
of th at has been questioned, most recentlv bv iMussack and Dappen l 201Cll . 



l201l[ ). In convectively unstable regions, the transport of heat is modelled by 
a (usually local) mixing-length formalism; Reynolds stresses are usually, but not 
always, ignored. Chemical species are presumed to be homogenized in convection 
zones; and elsewhere, aside from nuclear transmutation, chemical differentiation 
is purely by gravitational settling (and possibly radiative levitation) against mi- 
croscopic diffusion. Energy and material transport by the tachocline circulation 
is normally ignored, as is transport by acoustic and gravity waves generated 
by the turbulent convection. The total mass of the Sun is presumed to be con- 
served, and the "best" equation of state, opacity "formula" and nuclear reaction 
cross-sections are employed. 

It goes without saying that in computing hydrostatic structure it is prefer- 
able (some would say mandatory) that the coordinate singularity at the centre 
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of the star and the effective branch point at the radiative-convective interface 
(when a local mixing-length formalism is adopted) be correctly treated by the 
numerical integration, accurately enough for the resulting solar model to possess 
well defined seismic eigenfrequencies. Unfortunately, that is not always achieved, 
although I suspect that the numerical precision is greater than the observational 
accuracy. 



5. The Chemical Composition of the Sun 



Some time ago solar modelling {e.g. Christensen-Dalseaard et ai . 1996f ) had 
achieved fair agreement with the measured seismic structure of the Sun {e.g. 
Gough et Indeed, iBahcalll (j200lf ) considered that to be a "triumph 



for the theory of stellar evolution" , which he illustrated by plotting the relative 
devi ation of the solar sound speed from that of one of his standard solar models 
{e.g. iBahcall. Pinsonneault. and Basu . 2001 ). It is reproduced here as Figure 
[SI Yet might one not naively consider it to be an even greater triumph for 
helioseismology? After all, the seismological errors are much smaller than the 
modelling errors, as is more clearly evident in Figures [3] and 21 But I have already 
emphasized that the physics of seismology is very simple, whereas the physics 
of the structure of the Sun is not. The comparison can therefore be regarded in 
two opposite ways: i) one might marvel, as did Bahcall, that such complicated 
physics has been successfully reined to reproduce observation (more or less), or 
ii) one might consider it unremarkable that a mature theory, which has already 
survived a wide range of astronomical tests relating to other stars, and which 
depends on so many somewhat uncertain processes, can be adjusted to reproduce 
the measurable properties of the Sun. 1 think 1 side with Bahcall on this matter, 
because nobody has actually succeeded in reducing the discrepancy between the 
models and the Sun sufficiently to come within the much smaller seismological 
uncertainty. Moreover, the sit uation has been e xacer b ated by the results of recent 
spec troscopic re-ana l yses, by lAsplund et al. ( 2009[ ) , Caffau et al. (l2009l . I2OIII ) 



and Grevesse et al. 



( 201l[ ). of some of the chemical abundances in the solar 
atmosphere, those analyses now taking explicit account of spatial inhomogeneity 
in th e Sun's atmosphere caused bv convective motion as simulated bv Asplund 



et al. (l200nh based on a c ompressible hydrodynamical procedure described by 
Stein and Nordhmdl (Il998lfl. The earlv publications fAllende Prieto, Lambert , 



and As plund. I200ll.l2002l: lAsplund et all \2004 lAsplund. Grevesse. and Sauval 
2005al : lAsplund et all l2005bl lcf) suggested an enormous reduction in the total 
heavy-element abundance Z below previously accepted values: some 40 % or so. 



°The original analysis employed simulations in a model atmosphere with the previously ac- 
cepted chemical composition; from 2009 onwards, newer simulations were carried out with the 
Fe abundance inferred from the 2005 abundance determinations; the abundances of C, N, and 
O were not changed because their spectral lines are weak and have no significant impact on 
the radiative energy flux (R. Trampedach, personal communication, 2012). 
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If that represents the abundance in the radiative interior, it imphes, on average, 
a reduction in opacity of some 30%. By present-day standards, that figure is 
enormous. 



One of the best and most widely adopted solar models (combining all the gen- 
erally accepted physics for so-called "standard" theory into a model computed 
accurately enough for reliable computations of its oscillation eigenfrequencies to 
be possible) is the Model S discussed by IChristensen-Dalseaard et 'al. 
It was computed with a local mixing-length theory ignoring turbulent pressure, 
with mixing-length parameter [a] and initial hydrogen and heavy-element abun- 
dances [Xq = 0.7091, Zq = 0.0196 ] chosen to yield the correct present-day 
radius and what was believed to be the correct luminosity, and also a resulting 
present surface abundance ratio [Zg/Xs] of 0.0245, consistent with the spectro- 
scopic analysis of Grevesse and Noel^ (|l993 ). at an age tQ = 4.60 Gy (notwith- 
standing the different value quoted bv IChr istensen-Dalsga ard et all (|l996l ). and 
subsequently bv iDogan. Bonanno. and Ch ristensen-Dalsg aard ( 2010l )V The in- 
dividual surface abundances are Xg = 0.7373, Zg = 0.0181; the present surface 
helium abundance is = 0.2447, which agrees with observation within the 
limits of accuracy of the equation of state and helioseismological analyses of the 
depression of 71 caused by the second ionization of helium in the adiabatically 
stratified convection zon e (jGoughl. Il984bl: iKosovichev et al\. Il992l: Vorontsov. 
Baturin, and Pamvatnvkh. ll992l : lBasu and Antia . 1995 : Dappen . 2007 *). as I shall 
discuss briefly later. Figures |3] and |4] illustrate how well the seismic structure of 
the model corresponds to that of the Sun. 



The heavy-element abundance now favoured bv lGrevesse et al. (2011) is some- 
what greater than the original announcement, namely Zg = 0.0134 ± 0.0005; 
an independent and somewhat different analysis by ICaffau et al\ (|201ll ) that 
also takes convective inhomogeneity into account has yielded Zg = 0.0153 ± 
0.0011. Both of these values were estimated from abundances of the major 
opacity-producing elements except neon, whose abundance in the Sun cannot 
be i neasured spectroscopicallv. A recent compilation bv Lodders, Palme, and 
Gail ( 2OO9I I using a wider range of solar data, with an eve also on meteoritic 
abundances, yet using simple hydrostatics, has yielded the recommendation 
Zg/Xg = 0.0191. Accepting Yg = 0.249 from a calibration of model envelopes 
by iBasu and Antial ( 20041 ) using seismically accessible (asymptotic) integrals 
that are sensitive to the depression of 71 due to He 11 ionization, or = 0.224 



from a calibration of complete solar models (jHoudek and Goughl . l201lh . yields 
Zg ~ 0.0141 and Zg = 0.0145, respectively. (It is unclear whether calibrating 
complete solar models in this fashion is more reliable or less reliable that cal- 
ibrating only model envelopes.) These values are all substantially lower than 
that of Model S, on average by about 20%. What are the implications of these 
results? 
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5.1. The Abundances in Context 



Many of the discussions that have ensued have carried out solar-evolution com- 
putations and then compared, in some manner, implied seismic data, or seismic 
structure (principally sou nd speed), with the inferences from the Sun. They 



have been catalogued bv iBasu and Antial (|2008f ). However, seismic structure is 



merely a diagnostic of the issue, not the issue itself. It is instead much more 
lucid to the average physicist to be confronted with the direct implication of the 
newly reported abundances. So let us accept (most of) the assumptions of solar 
modelling, as is usually done in the discussions of the problem, and ask what 
they really imply. I start with the most obvious: that the Sun approached the 
main sequence fully mixed, having just come down the Hayashi track. There- 
fore the heavy-element abundance [Z(r, t©)] in the radiative interior today is 
(almost) the same as it is in the photosphere, aside from a relative excess of 
3 % or so due to differential gravitational settling. Heavy elements influence the 
structure of the Sun primarily via their effect on the opacity [k] , which controls 
the relation between luminosity [L] and temperature [T] through the equation 
of radiative transfer: L{r,t) ~ — (iQnacr^T^ /3kp^ dT/dr. However, T is not a 
seismic variable, and therefore cannot be measured directly. It is related to the 
seismic variables and p by the equation of state: T = T {c^ , p;Y, Z) - I set 
aside, for the moment, consideration of the relative abundances of the heavy 
elements. The contribution of Z to the equation of state is only about 0.5 %, so 
for the purposes of this discussion its uncertainty can safely be ignored. There 
remains only the helium abundance [Y] which is to be estimated from the theory 
of stellar evolution. It is important to appreciate that this is the only stage in 
the argument where the details of the evolution theory come into serious play. 
However, it must be realized that knowledge of the functional form of Y{r) 
today is crucial for assessing the stratification of the Sun's radiative interior. 
Therefore, some attention must be paid to how it is determined, appreciating 
the assumptions to which the outcome is sensitive. 

As the Sun evolves on the main sequence, hydrogen is transmuted into helium 
in the core, adding a spatially and temporally varying component df^Y to the 
helium abundance: Y{r,t) = [1 + sy(r, t)] Yq + S^Y, where syFo is the change 
in Y produced by g ravita tional settling moderated by diffusion. It is known 
{e.g. iGoughl . Il983d . [ 1990( ) that the time dependence of the total luminosity 
Ls{t) L{R,t) of the solar model, calibrated to satisfy I/s(io) = Lq, is not 
very sensitive to the details of the theory (such as the choice of Zq, or whether 
or not there has been some small degree of material mixing in the core - we 
know from seismology that the core cannot have been homogenized, a conclusion 
which is consistent with HR diagrams of solar-like stars), so the total amount 
of hehum that has been produced by today [/ JeFdr] is proportional to / Lgdt, 
which is essentially known (provided Iq is known). Therefore it is adequate for 
the current discussion to accept the function S^Y from any standard model. Also, 
the function sy(r, i) is only very weakly dependent on Yq, and may also safely 
be taken as given by the model. Likewise, one can take Z = [1 + sz{r,t)] Zq 
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with sz given by the model (it is broadly similar to the function sy). Whence 
T(r, t) =T {c^,P', (1 + sy) Fo + ^6^7 (1 + Sz) ^o) , where and p have the values 
determined seismologically from the Sun. Only the initial abundance Yq now 
remains unknown. How do we determine it? 



I recommend calibrating Yq by accepting the nuclear reaction rates - after 
all, all the pertinent nuclear cross-sections have been carefully assessed experi- 
mentally during the investigations of the neutrino problem, save the controlling 
p~p reaction cross-section which is determined only theoretically (however, it is 
the simplest of all nuclear reactions). There remains also an i ssue concer ning the 
screening of energetic particles (e.g. iMussack and Dappenl . l2010l . l201lh . These 
matters are probably, for the purposes of this discussion, minor. Therefore the 
total rate of thermal energy production, hence the luminosity, is determined in 
terms of and Yq. Granted that Lq and t© are known quite well, Yq can thus 
be calibrated to determine T(r) todajH, which can then be substituted into the 
equation of radiative transfer to evaluate K{r). That procedure can either be 
carried out explicitly from a seismological inversion to determine the seismic 
stratification, or imp l icitly within the inversion itself, such as in the manner 
advocated by lEUiottl (Il995l) . The outcome is illustrated in Figure [BJ in which 
the continuous curve is the relative difference {kq — ks)/ks between the Sun's 
opacity [kq] determined in the manner that I have just described and that of 
Christensen-Dalsgaard's Model S. Not surprisingly, kq is very similar to the 
opacity in Model S. The issue posed by Asplund et al. and Caffau et al. is 



^Impl icatio ns from earlier estimates of Y to determin e T{r) have been discussed by lEUiot j 
1 II995D and iTripathv and Christensen-Dalsgaardl lll998f) , the first under the assumption that 
Y(r) differs from that in a reference solar model by just a constant, the second that it can 
be obtained simply by scaling the reference-model value by a constant factor, both of which 
are inaccurate in the core, although Elliott notes that the opacity perturbations produce a 
predominantly local response, so that these analyses should provide a fair estimate in the 
radiative envelope. The Triga^i^^hnsten^er^Dalsgaa^^c^in^^a^ise^^^Trip Basu, 
and Christensen-Dalsgaard to obtain the opacity difference from Model S by RLS 

frequency fitting (J. Christensen-Dalsgaard, personal communication, 2012) assuming that that 
difference can be expressed as a function of T alone, yielding a superficially similar functional 
form to the continuo us curve in Figure 6, but with a magnitude about 50 per cent greater. 
iBahcall et al\ | |2005D have estimated the opacity difference by adjusting k by hand in solar 
models. Their preferred model had a const ant 11 % augm entation over the OPAL valu es using 
the most recent abundance determinationsj^^glun^^^i. , 2000, 200 4: lAsplundl. 120051: AUende 
Prieto, Lambert, and Asplund' T200iri200^ !Asphm^Tersonal communication with Bahcall et 
al., 2004) in the radiative envelope down to T = 5 X lO^K, beneath which the augmentation 
was smoothly reduced to zero (probably by a half Lorentzian function with half- width at half 
maximum of 2 X lO^K), and has a seismic structure as close to that of the Sun as does Model 
S. Regarding that model as a proxy Sun, one would expect the relative opacity difference 
be tween it and a model with the unm odified abundances to be comparable with the inference 
by IChris tcnscn-Dalsgaard_ef_ajjj200g[^^JIlieJ^ Korz ennik 
and U lrich (J.989) had earlier estimated opacity errors by RLS (L2 norm) data fitting, and lSaiol 
1I1992D by Li data fitting, each by scaling k by a function of T and ignoring the dependence of 
the relation between T and the seismic variables on chemical composition; they expressed their 
results as deviations from different reference model s, so they cannot easily be compared with 
those presented in Figure 6. An estimate by OLA llTakata and Goughl . 120011 ) of the absolute 
structure of the Sun, including opacity, using the procedure for de termining Y described in 
the te xt (together with tachocline homogenization as c alibrated bv[Elliottj Gough, and Sekiil 
lll99Sl) 1 is presented bv lGough and Scherreil ||2002|) and lGoughl ||2004 l20odV 
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Figure 6. The dot-dashed curve is the estimate by IChristensen-Da lsgaard et al\ ll2009f) 
of the relative opaeity augmentation require d when adopting the |AsDhmdril2005l) chemical 
composition in OPAL opacity computations l llglesias and Roger s!. IT99^ to produce a solar 
model with the same sound speed as in Model S of lChristen scn-Dalsgaard et a l\ ||1996[). The 
dashe d curve represents the augmentation, compute d by Ic hristens en-Dalsgaard and HoudeS 
1I2OIOD . required of OPAL opacities with the revised lAspIund et al .1 II2009D abundances. The 
triple-dot-dashed c urve is a comparable difference between two solar models computed by 
iBahcall et al\ | |2005D with OPAL opacities, the first with the latest abundances by Asplund et 
al. at the time, the second with abundances that have been artificially chosen to yield what 
one might regard as a seismically acceptable proxy Sun. For comparison, the dotted curve is 
a linearized estimate of the same quantity were the relative abundances of the heavy elements 
to have beer^nrese^ed^isin^^gagit^^firiva^ve^^bteine^^o^ and 
Pinsonneault 1 II992D . The continuous curve is the relative difference (kq — Kg)/Ks inferred by 
iGoughl ^Mi) between the Sun's opacity [kq\ and that of Model S. 



therefore simply an abundance-opacity problem: how can ^©(r) be reconciled 
with their abundance measurements? 



5.2. Suggestions for Reconciliation 



An obvious naive suggestion is that opacity calculations are in error by just 
the appropriate factor required to compensate for the proposed revision in Z, 
leaving the functional dependence on the other variables unscathed. Given the 
testing that has been undertaken by Iglesias and Rogers, that does not seem 
so very likely. It is also unlikely that energy transport is opposed by some 
other mechanism, such as gravity-wave transport, because that would require 
a completely different physical process to mimic the functional form of the effect 
of radiatively induced atomic transitions. Of course, it could be that the new 
abundance determinations are in error, despite the extra care that has been 
taken; we ha ve seen an increase since the earlv work by AUende Prieto, Lambert, 
and Asplund lli'oOll . |2002[ ) and I Asplund et al\ {20^ : Asplund, Grevesse, and 
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Sauval (|2005ah : lAsDlund et all (|2005bl R). so a further increase might not be too 
surprising. However, the independent investigation by Caffau et al. (12009 .1 2011 ) 
has yielded results no t very dissimilar from those of lAsplund et am 20091 ) and 
Grevesse et al\ (|201l[ ). which adds to the crcdibihty of the new abundances. Of 
course, another possibihty is that early in its post-Hayashi days the Su n was 
contaminated by metal-poor material, as iGuzik. Watson, and Co^ (|2005l ) have 
suggested. If that were the case, it may just be possible to detect seismologically 
a relic compositional discontinuity at the location of the base of the zero-age 
convection zone, provided that it has not been destroyed. The discontinuity is 
at least stable to double-diffusive convection, so it would not be destroyed spon- 
taneously. Yet a fourth possibility is that the relative abundances of the heavy 
elements are different from what is gene rally believed, p articularly that of neon, 
as has been discussed, for example, bv lOrake and Test a (2005^), leaving open 
the possibility that Z is greater than the values suggested by Asplund et al. and 
Caffau et al. based on the assumption that the relation of the neon abundance to 
the abundances of the other opacity-producing elements is preserved, and that 
therefore k is perhaps proportionately greater too. That too is regarded by many 
spectral analysts as being intrinsically unlikely, and in any case neon alone does 
not mimic the heavy-element mixture adequately to restore the solar models to 
their pre- Asplund state. 



The apparent opacity discr e pancy has evolved with the abundance revisions 

(^2009^ estimated the amo unt S k by wh i ch the 



Christensen-Dalsgaard et al. 



opacity according to OPAL (jiglesias and Roeersl . 119961 1 with lAsplundl (|2005l l 
chemical composition would need to be augmented in order to produce an 
otherwise standard solar model with the same sound speed as in Model S, 
and hence about the same sound speed as in the Sun. It is depicted as the 
dot-dashed curve in Figure El where it is plotted against log(T) along the 
thermodynamic p-T path of one of the models. One might have expected the 
outcome to have been given approximately by {d In n/d In Z)p T,x S In Z, in 
which S In Z In Zs — In Z\sp ~ In .^s.s ~ In -^s.Asp = 0-370, the subscript s 
denoting surface value and the subscripts S and Asp denoting Model S ( Grevesse 



and Noels. 1199,'^ and lAsplundl (l200fJ ) values respectively, and where the opacity 
derivative is evaluated at constant relative heavy-element abundances. However, 
it is evident from Figure IH] that that is not the case. According to J0rgen 
Christensen-Dalsgaard (personal communication, 2011) the quite substantial dif- 
ference arises because the relative heavy-element abundances in the two models 
differ. Unfortunately, that renders back-of-the-envelop e estimates suspect. The 
later abundance determinations bv I Asplund et al. ( 2009t ) have brought the OPAL 
opacities closer to th ose of the Sun, a s the dashed curve in Figure 6 depicts. 
It is interest ing that Antia and Basul (l201ll ). using the chemical composition 
proposed bv lCaffau et al\ (12011), whose relat i ve abu ndanc es do differ from the 
earlier values given by iGrevesse and Sauvall (Il998l) . and IGrevesse and Noeld 



(|l993r ). and presumably also the more modern values of IGrevesse et ^ daonl) 



report that they have obtained solar models t hat are seismically almost as good 
as models constructed with the abundances of IGrevesse and Sauvall (1 19981 ). and 
they found even closer agreement with the Sun if the assumed abundance of neon 
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were artificially enhanced by a factor That essentially reduces the problem 
of reproducing the Sun's seismic stratification theoretically to the state in which 
it was prior to Asplund's original announcement. It must be realized, however, 
that that does not close the matter, because merely reproducing previous partial 
results does not necessarily prove the veracity of the physics behind the new 
models. Moreover, the seismic structures of the models, new and old, deviate 
from that of the Sun by many standard errors, as is evinced by Figure [31 



5.3. Possible Flaws in the Argument 



The foregoing discussion is predicated on the presumption that the sole signifi- 
cantly discrepant ingredient in the solar modelling is opacity. That need not be 
the case; modification of any non-seismic variable might, at least in principle, 
bring a theoretical model into line with observation. Evidently, a direct change in 
the relation between the seismic variables pressure and density, for example, via a 
modification to the abundances of the abundant chemical elements H and He, can 
alter the sound speed: Gough and Kosovichev ( 1988 . 19901 ) found that the strat- 
ification of the energy-generating core can be reproduced by a slight smoothing 
of the abundance profile that had been produced by nuclear transmutation, sug- 
gesting that a slight degree of mixing has taken place, possibly by disturbances 
that have been shear-generated in a manner analogous to clear-air turbulence 
in the Earth's atmosphere. Alternatively, there could be an additional source of 
energy transport without mixing, such as by accreted weakly interacting massive 
particles (wimps) , whose presence in the Sun was postulated originally in order 
to try to account for the observed low neutrino flux without neutrino transitions 
(jSpergel and Press! . llQSST ) ; wimps modify the temperature distribution, and con- 
sequently the distribution of helium produced by nuclear transmutation, leading 
to a sound-speed modification whose general functional form is broadly similar 
to what would be required (jGilliland et alV fl98ftl . although a serious attempt 



to reproduce the seismic structure of the radiative interior appears not to have 
been made (see also Faulkner. Gough. and Vahia ^ 19861: Da ppen, Gilliland, and 
Christensen-Dalsgaard, 1986 : Christensen-Dalsgaardl . 1992 ). 



It is worth commenting that a reduction in Zg at fixed relative abundances 
of only about the 20 % below that of Model S implies, via the theory of solar 
evolution, a reduction of about 0.02 in Kj- This appears to me to be not dissimilar 
to the margins of the true uncertainty in the helioseismological determinations of 
the helium abundance of the convection zone - although not within the precision 
quot ed by those who have attempted to determine by only a single procedure 
(e.^. |Basu and Antial. Il995l. l2004l: iBasul. Il998l: iRichard et ali Il998l: Pi Mauro 
et ai. 120021 ). even granted that systematic errors had been recognized. 
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5.4. Seismological Investigation 



Chemical abundances can in principle be measured directly by seismology by 
analysing the non-ideal properties of the solar plasma. Specifically, 71 is de- 
pressed in regions of partial ionization, by an amount that is almost proportional 
to the abundance of the ionizing clement. IPappen and Gough ( 1984f ) proposed 
measuring the helium abundance by calibrating the ionization-induced variation 
of a thermodynamic function [Q] expressible in terms of 71 and its derivatives and 
which is easily accessible to seismological probing in the adiabatically stratified 



layers of the convection zone (jGoughl . Il984bl . A procedure for so doing was 
developed bv lPappen. Gough. and ThompsonI ( 1988f ). who found Yg to be lower 
than tvpical values obtained from calibrating evolved solar models (Dappen 



et al., [I991D, which triggered IChristcnscn-Dalsgaard. Proffitt. and Thompson 
(I1993I I ;o investigate the seismological implications of gravitational settling of 



heavy elements, with notable success. Other forms of calibration, often based 



more directly on 71, have been pursued {e.g. Dappe n et aLl . ll991l : lBasu and Antia 



19951 . l2004l:rBasul.ll998l: [Richard et al\.\l998: Pi Mauro et a/.l. I2OO2I: Houdek and 



Gough, 120071 ). None has been completely satisfactory, not least because there 
are serious uncertainties in the equation of state. There have been attempts to 
estimate Z by calibrating solar models against low-degree frequency-separation 
ratios assuming tQ to be known {e.g. Chaplin et al . 20071 ): that is equivalent to 
calibrating Y , of course. 



A direct seismological estimate of Z in the convection zone, and thus a direct 
estimate of Zg, via the ionization-induced depressions in 71 or in the variation of 
the seismologically more accessible function of 71 and its derivatives, is very 
delicate. It would probably be necessary to average the individual variations 
in the different ionization zones, so the result would depend to some degree 
on the assumed values of the relative abundances ( Mussack and Goughl . [2009h . 
Nevertheless, if successful, the accuracy of the result could probably be assessed 
from the seismic frequency uncertainties, because the relation between the mag- 
nitude of the 71 depression and the abundance of any ionizing element producing 
it, at least at the level required to judge the abundance-opacity problem, is 
fairly robust. The seismologically apparently more straightforward procedure of 
measuring a spatially averaged absolute value of and relating that to the 
sligh tlv greater value expected in the absence of heavv elements, as Antia and 
Basu ( 20061 ) have tried, is more uncertain because the change in 6 due to heavy- 
element ionization that is sought appears to be less than the uncertainty in the 
equation of state arising, for example, from the commonly neglected, at least 
in the chemical picture ( Dappenl . 12004 . 2007 ). finite volume occupied by bound 
species {e.g. Baturin et al. . 2000l ). 
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5.5. The Equation of State 



My statements about the uncertainty in Yg and the equation of state call for some 
justification. It goes without saying that chemical composition is not a seismic 
variable, not even the abundances X and Y of the abundant elements H and 
He. In order to determine say, it is evidently necessary to use an equation of 
state, which relates the seismic variables to composition via another, necessarily 
non-seismic, variable, such as T or s, and subsequently, at some point, to confine 
attention to the adiabatically stratified region of the convection zone where that 
variable can be eliminated from the spatial variation of the seismic variables. 
Once that has been accomplished one can learn about chemical composition 
through the effect of ionization on 71 . This has been approached either directly 
by measuring some appropriate property of 71 itself, or by working with some 
thermodynamic function 8 of it. The reliability of the outcome necessarily rests 
on the reliability of the equation of state, which is very difficult to assess. That 
should be obvious because there is no redundancy in the dependence of the 
seismic modes on the s eismic variables. Therefore, as has been alluded to in the 
past {e.g. Gough . 2004 ). perhaps too obliquely, any intrinsic error in the equation 
of state cannot be assessed by seismology alone. To make progress, non-seismic 
information must be incorporated. In practice, that information comes from 
some prior appreciation of the reliability of some aspects of the equation of state, 
coupled with the additional non-seismic constraint that the convection zone is 
chemically homogeneous (sometimes augmented with the constraint that deep 
down the stratification is adiabatic). Suppose, for example, one is estimating the 
deviation of the Sun from a reference solar model, represented by (Sln7i and Slnu, 
where u{r) is a complementary seismic variable. From seismology one can relate 
averages (51n7i +(51nu of those deviations to data combinations d, where the single 
and double overbars indicate simply that the averaging kernels are different. 
Then, one can write the first average as the sum of {dlnji / d\nu)Y 6\nu, which 
can be incorporated into the second average, of {dlnji /dY)u SY, where 5Y is the 
difference in Y between the Sun and the reference model, and of a component 
i5intln7i resulting from intrinsic error in the equation of state. One would like 
to be able to distinguish between those components. However, that is strictly 
impossible because they have the same averaging kernel. Therefore one cannot 
unambiguously determine SY from the data [d] . What has been attempted in the 
past is to design kernels such that those defining the overbar and the double over- 
bar are both in some sense small, and then to estimate SY by neglecting (Sintln7i 
and Slmi, a procedure which is evidently not strictly valid; after that, (5intln7i can 
be computed with presumably different and well lo calized averaging kernels. The 
results of such a procedure ha ve be en presented bylBasu. Dappen. and Navfonov 
(|l999l ). lDi Mauro et all tOOi ). and loappenl (|2004h . using both OPAL and MHD, 
the two most popular, and probably the best, equations of state available today. 
Not only do their inferred intrinsic errors from the two equations differ, as they 
surely must, but so too do the estimates of the helium abundance [Yg] in the Sun's 
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convection zon43- It must be appreciated that any seismological inference about 
'5intln7i such as this is susceptible to potential errors in the inferred values of Yg, 
which are difficult to appreciate because the precise manner in which SY was 
obtained is unclear. However, well outside the hydrogen and helium ionization 
zones dliryi/dY is small, so errors in 6Y can hardly contaminate (5intln7i there. 
This is a robust property of any realistic equation of state, because under normal 
stellar conditions 71 is essentially independent of Y where hydrogen and helium 
are both fully ionized, errors in the small dependency that does remain (arising 
from processes such as the influence of the electron density on the ionization 
of heavy elements) being characterized by the quantity <5intln7i itself. However, 
the situation is different within the ionization zones, where in general (Sintln7i 
and (91n7i /dY) SY contribute comparably, and cannot be entirely recognized 
apart. I must add, however, that Slnji and (dhi'ji/dY) are normally functionally 
different, so 5intln7i cannot vanish everywhere. 



An alternative, more transparent, approach to an attempt at separation could 
be to adopt the attitude that one's chosen equation of state is as good as it can 
be, at least in the adiabatically stratified region where Shryi has been inferred. 
Therefore one could choose for each equation of state the value of 6Y that 
minimizes the integral with respect to acoustic radius of ((5intln7i)^ over that 
region - I choose acoustic radius because that is the natural seismic frequency- 
controlling independent_vBjiaMB. outcome, using the inferences of (Sintln7i 
and Fs bv iDi Mauro et all (|2002l l. is Y^ = 0.2430 for OPAL and Y^ = 0.2280 
for MHD; I do not quote a precision because it is hardly material to estimating 
accuracy in this context. The range of all these values, roughly 0.02, must surely 
be regarded as a lower bound to the range of values within which the solar Ys 
is likely to lie, because the discrepancies are systematic; that value is some 40- 
100 times gr eater than tvpical errors quoted bv. for example. iBasul (1998) or Di 
Mauro et al. (I2OO2I) . based on the precision of the specific procedures that were 
carried out. Based solely on the calibrations described in this paragraph, one 
might optimistically conclude that Yg = 0.24 ± 0.01, while recognizing that the 
accuracy might well have been overestimated. Consequent estimates of (Sintln7i 
are illustrated in Figure [T] 



■'Di Mauro et al. quote Ys = 0.2539 ± 0.0005 when OPAL is used, Ys = 0.2457 ± 0.0005 
when MHD is used, the errors being merely formal, representing a precision error that takes 
no account of the error in the relation between Ys and 71 in the reference model. With those 
values the magnitude of the inferred 5intln7i was found to be the greater for the MHD equation 
of state beneath the helium ionization zones, and the lesser above r/R 0.97 in the Hell 
ionization zone. That is perhaps not surprising because MHD is possibly better at taking into 
account the complicated chemistry that dominates higher up in the solar envelope, whereas 
the virial expansion used for OPAL is perhaps more reliable where such complications make 
only a minor non-ideal contribution. The difference between OPAL and MHD must surely offer 
some estimate of the total uncertainty. 
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Figure 7. Inferences concerning the intrinsic error in the OPAL (asterisks) and MHD (trian- 
gles) equations of state, based on the assumption that the hehum abundance, in each case, is 
such as to minimize its L2 norm with respect to acoustic radius over the region of the convection 
zone that is stratified adiabatically to better than 1 part in 10"^, computed from the analysis of 
iDi Mauro et al\ ll2002t) . Because the helium abundance is uncertain, the symbols should not be 
meant to represent reliable estimates of the actual intrinsic errors 5i,itln7i. Formal error bars 
are not drawn because they are enormously smaller than the true uncertainties, and might 
therefore confuse. The dotted curve is 0.02(91n7i/0y)u, computed with the OPAL equation 
of state. 



5.6. Other Matters to Consider 



There are uncertainties in other quantities used in estimations of the chemical 
composition that still need to be addressed. Amongst them are the current 
luminosity [Lq] and the age [to], a matter to which I have already alluded. The 
luminosity is inferred from the total solar irradiance, and is currently obtained 
by assu ming that the radiant output is spherically symmetric (e.g. Willson and 
Hudson. Il988[ ). That assumption is n ot strictly correct, as models of solar- 
cycle irradianc e variation indicate (e.g. Taylor et ai . 1998 : Foukal et al. . 20061 : 



Frohlichl . I2OIII ): accounting for the asphericity augments the estimate of Lq by 
some 1.5% or so. However, there has been a recent downward revision in the 
measured value of the i rradiance from an end-to-end instrumental recalibration 
( Kopp and Leanl . l201l[ ). which is partially compensating. The age [t©] is more 
uncertain. It is likely to be about the same or only slightly greater than the age 
of th e oldest meteorites, which seem to lie between 4.563 and 4.576 Gv f Amelin 



et aL. l2002l : |jacobsen et ad |2008| . l2009l : iBouvier and Wadhwal . l2010f ). Attempts 
to calibrate solar models seismically have yielded values equal to or slightly 
lower than the age ts adopted for Christensen-Dalsgaard's model S, namely 
4.60 Gy. Those calibrations in which a va lue of Zg is assumed to be that of 
Model S yield a valu e rather lower than ^c; (Dziembo wski e t all 19991: Bonanno , 
Schlattl, and Paterno, 20021 : Dogan. Bonanno. and Ch ristcn sen-Dalsgaard . 2010( ). 
the most recent, based on BiSON data ( Chaplin et all i2Q07t ) being 4.57 Gy; 
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a simultaneous calibration of Zg and t© by iHoudek and Gough ( 2011 ). also 



with BiSON data (jBasu et all 120071 ). has yielded precisely ts- Interestingly, 



the heavy-element abundance in the latter investigation was calibrated to be 
Zs = 0.0142, es sentiallv the sam e as th e recent value preferred bv Loddcrs, 
Palnie. and Gail (|2009f). and Iving between the values preferred bv Asplund 



et al. ( 2009f ) and lCaffau et al. ( 2011 ). However, the model was calculated with the 
relative heavy-element abundances of Grevesse and Noel^ ( 19931 ) (neon not being 
artificially enhanced); therefore the opacity is too low and the seismic structure 
cannot be correctQ. Whether that has yielded a superior or an inferior estimate 
of <0 is unclear (jGoueh . 2012a ). The calibration relies partially on estimating 
1^ seismologically from the oscillatory component of the low-degree frequency 
distribution caused by the acoustic glitch associated with the depression of 71 
in the helium ionization zones (jHoudek and Goughl . 120071 ) . To accomplish that 
estimate it was assumed that hydrostatic support is solely a balance between 
pressure gradient and gravity, as is usual: magnetic stresses and the centrifugal 
force were ignored. (It was also assumed that the Sun is spherically symmetric. 



as in other calibrations, which we know is not strictly correct (e.g. iGoode et al 



2OO2I : iKomm. Howe, and Hil3 . I2OO6I : lEmilio et ail . l2007l : iFivian et all I2OO8 



Kuhn et al. . 2012f ): by how much the structural asymmetry contaminates the 
assumed relation between the degree-dependence of the seismic frequencies and 
the locations of their lower turning points, of crucial importance to the structural 
calibrations, is yet to be ascertained.) 

The validity of ignoring magnetic stress has been implicitly questioned by 
Basu and Mandel (I2OO4I ) and IVerner. Chaplin, and ElsworthI (|2006f ) , who re- 
ported a 10% solar-cycle variation in the glitch signal. That variation is huge, 
and, if correct, cannot possibly result from a temporal abundance variation. It 
has been pointed out that were the glitch variation to be magnetic, an intensity 
variation of some ten or so Tesla in the second He ionization zone would be im- 
plied ( Gough . 2006f ). Moreover, it would introduce an error in the two-parameter 
age calibration of similar relative magnitude ([Houdek and Goughl . 120111 ) . How- 
Christensen-Dalsgaard et all { 2011 ) have recently failed to detect a change 



of such magnitude, although they admit that further analysis is necessary in 
order to be sure. Meanwhile, iBasu et al\ (|201G| ) maintain that the change does 
occur. The whole issue is evidently very important, and we would like to know 
the answer. It may be some time before we do. 



^Moreover, the initial helium abundance [Yq] of the calibrated model is about 0.250, which is 
dangerously close to the amount [Yp] believed to have been created by Big-Bang nucleosynthe- 
sis, whose estimated value has been climbing over the last two decades (Stcigman, 2 0071 ): the 
latest estimate is Yp = 0.2478 ± 0.0006 (G. Stcigman and M. Pettini, personal communication, 
2011). Subsequent contamination of the interstellar medium by supernovae exacerbates the 
situation. 
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6. What of Importance is There Left to be Learnt? 



Helioscismology was amazingly successful in the early days. It demonstrated that 
the resolution of the solar neutrino problem was not to be found in adjustments 
to parameters of a standard solar model, and later, as precision increased, that 
the resolution must be sought in nuclear or particle physics. As we all now know, 
the matter is resolved by neutrino transitions. Helioseismology has also provided 
estimates of the solar helium abundance in the convection zone, perhaps not 
as accurately as we would like, yet which, granted gravitational settling, are 
at least not seriously at variance with Big-Bang nucleosynthesis. Coupled with 
that is a precise measure of the location of the base of the convection zon e 
( Christensen-Dalseaard. Gough. and Thompsonl . 1991 : Basu and Antia . 1997 ). 



which has been us ed extensively as a si r nple d i agnostic for calibrating theoretical 
solar mod els (e.g. Turck -Chicze et all 1993 : Bahcall. Pinsonneault. and Basul . 



[2001. :. Guzik. Watson, and C ox. 2005), and is important for defining the boundary 
conditions in numerical simulations of the convection zone. In addition, we know 
that the quadrupole moment [J2] of the ext erior gravitational e quiptentials is 



about 2.2 X lO^'^ fe.o. lSchou et aLl . Il998. : , Antia. Chitre. and Gouah . .200S ). con- 



tributing to the precession of planetary orbits, particularly that of Mercury, by 
an amount that is compatible with General Relativity. The latitudinal variation 
of the angular velocity [fi] in the photosphere persists, approximately, to the 
base of the convection zone, beneath which is a thin interface, the tachocline, 
and a rigidly rotating interior (admittedly with some doubt about the rotation 
of the energy-generating core) . 

We seem to have answered the major outstanding questions accessible to 
seismology that interest most of the scientific community. So is global helioseis- 
mology finished? Is helioseismology finished? To be sure, there is the unfinished 
business I talked about in the previous section, but it may require only minor 
details to sort that out. It is now often suggested by commentators that the days 
of true excitement in the subject, and of important discovery, are over. 

I dissent vehemently from that view. Firstly, the issues that I discussed in 
the previous section are related to the microphysics of opacity and the equation 
of state, which are important to plasma physics and may end up being very 
important also to the study of stars other than the Sun. More detailed seis- 
mological investigation may be able to unveil, or at least provide clues to, some 
fundamental unresolved issues in plasma physics, such as, in the chemical picture, 
how the energy difference between the continuum of "unbound" electrons and 
the strongly bound states is determined. This is related to the whole matter of 
electron screening of charged species, a reliable consistent quantum-mechanical 
study of which is still wanting, a nd has implications also for thermonuclear 



reaction rates {e.g. iDappenl . Il998f ). But more obvious is the whole matter of 



the internal macroscopic dynamics of the Sun: large- and medium-scale angular- 
momentum transport; material redistribution by meridional flow in the radiative 
zone; the pattern of the large-scale meridional flow in the convection zone; 
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augmentation, distortion, and decay of the magnetic field and its back-reaction 
on the flow; the formation, evolution, and final decay of sunspots and other 
forms of activity; and how all these conspire to drive and control the solar cycle, 
and, perhaps most prominently, modulate the solar outputs of particles and 
electromagnetic radiation that influence the climate on Earth. Many of these 
are the object of ongoing local hclioseismological techniques - time-distance and 
seismic holography, and ring analysis. I have not discussed those techniques here. 
Ring analysis is relatively straightforward, and has given us views of horizontal 
flow not far beneath the photosphere. Telechronoseismology has given us a view 
of a sunspot, complete with the flow around it: deep divergent horizontal flow 
carrying away the excess heat rising around the obstructing spot and an as- 
sociated convergent subohotosDheric counterflow above fZhao, Kosovichev, and 
Duvall, 2 001; Zhao and Kpsovichev,, 2003; Kosovichev. ,2009i: Zhao. Kosovichev, 
and Sekii, 2010[ ). Di vergent (Evershed ) motion has now been detected even closer 
to the photosphere ( Zhao et aZ.I . l20li[ ) - a critical feature that was missing from 
the original scismological analysis, and which caused some onlookers to harbour 
grave doubts. The qualitative picture, which has made its way onto the well- 
known SOI/MDI coffee mug, is now not only physically plausible, but also in 
accord with photosphe ric observations. A dynamical picture of the Eversh ed flow 
is now emerging {e.g. IWeiss et all. l2004l: Kitiashvih era/]. l2009l . l2010h . More 
recent helioseismic observations ( Kosovichev and Duvalll . 12011 ) have revealed 
the cessation of the deeper subphotospheric convergent flows as sunspots decay, 
flows which one might presume had previously held the spots intact. However, 
different scismological analvscs arc not vet all entirelv consistent (e.g. Gizon 
et a/.~l2009; Hindman, Haber, and Toomrel . bOP'of T so the picture is not wholly 
secure. Also, more extensive observations will be needed for understanding the 
entire life-cycle of the spots. 



Another matter of importance for understanding the global internal dynamics 
is the structure of the tachocline and its associated meridional flow, and how 
that influences the convection above. It is extremely difficult to detect deep 
meridional flow seismologically, because the frequency perturbations are small; 
however, so me progr ess ap pears to be possib le with the use of eigenfun ction 
distortions ( Schou. Woodard. and Birch . ,2009i : Gough and Hindman , 2010( ) , and 
lAntia. Chitre. and Gough ( 2012f ) are brving to advance simplified seismolog- 
ically constrained dynamical arguments. If the meridional tachocline flow is 
downwelling near the Equa tor and the Poles , and u pwel ling near the latitude 



of zer o tachocline shear, as ISpiegel and ZahnI (jl992l ) and iGough and Mclntvre 



(jl998f ) have argued, how does that flow connect to the general meridional flow 
in the convection zone? Or is that flow so slow that it is simply swept aside 
unnoticed by an independent Reynolds-stress-driven circulation? Does the up- 
welling dredge up a primordial magnetic field from beneath, and if so, does the 
shearing of that field in the tachocline react back on the rotation to produce a 
mid-latitude shear-free region? Does the dre dged field stoke a (non)dynamo in 
the convection zone (^cf. iBvington et at. . 2012)? Does the equatorward tachocline 
flow in the polar regions and the poleward flow in the equatorial regions advect 
the periphery of the predominantly dipole remnant field in the radiative interior 
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to align with the upwehiiig flow- convergence zone (Gough, l2012cf ). or do the 
torques from the rotational shear dominate the dynamics, as Toby Wood and 
Michael Mclntyrc (in preparation) quite plausibly presume? Is such a putative 
inclined dipole responsible for the active longitudes? - and the emergence of 
sunspots? Most previous dynamical studies have presumed the tachocline flow 
to be csscntial lv steadv: the evidence presented recentlv bv Antia. Chitrc, and 
Gough (j20I2,) suggests that it varies with the solar cycle. So these issues will need 
to be readdressed in a new light. The rise of su nspots through the convection is 
now b ecoming accessible to seismic observation (|llonidis. Zhao, and Kosovichevl . 
201l[ ). There arc also questions related to structural changes associated with 



the cycle, and how they arc related to the variations in the Sun's outputs. And 
what of the rotation of the solar core, and the associated circulation, if 51 differs 
substantially from the angular velocity of the surrounding envelope? That is 
likely to influence our ideas about other stars. 



It must be appreciated that the Sun is an important benchmark for the whole 
of stellar physics, and helioseismology for asteroseismology. Finding answers to 
many of these questions is likely to be assisted by new helioseismological findings, 
although seismology alone will not be enough. Answers will not come easily - 
the easy questions were answered quickly (although not necessarily easily) in the 
early days. But the rewards from answering the new questions - designing subtle, 
hardly detectable, seismic diagnostics and developing techniques to analyse them 
- are potentially great. Their pursuit is to a large degree the task of the new 
generation of young seismologists. 
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